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The glass transition, extensively studied in dense fluids, polymers, or colloids, corresponds to a 
dramatic evolution of equilibrium transport coefficients upon a modest change of control parameter, 
like temperature or pressure. A similar phenomenology is found in many systems evolving far 
from equilibrium, such as driven granular media, active and living matter. While many theories 
compete to describe the glass transition at thermal equilibrium, very little is understood far from 
equilibrium. Here, we solve the dynamics of a specific, yet representative, class of glass models in 
the presence of nonthermal driving forces and energy dissipation, and show that a dynamic arrest 
can take place in these nonequilibrium conditions. While the location of the transition depends on 
the specifics of the driving mechanisms, important features of the glassy dynamics are insensitive 
to details, suggesting that an 'effective' thermal dynamics generically emerges at long time scales in 
nonequilibrium systems close to dynamic arrest. 



Our goal is to study theoretically whether a glass 
transition can occur in a many-body interacting system 
driven by fluctuations which do not have a thermal ori- 
gin. We motivate this very general question using recent 
experimental findings, which we divide into two broad 
classes. First, in the presence of mechanical forcing, gran- 
ular media display a form of dynamical slowing down 
which shares important similarities with the dynamics 
of molecular liquids near the glass transition. Forces 
driving the dynamics of the grains can for instance be 
a global oscillatory shear high-frequency vertical vi- 
brations J! or an air flow through the granular as- 
sembly [4]. The injected mechanical energy produces er- 
ratic grain motions, while energy is dissipated through 
collisions. Dense granular fluids driven through dynamic 
arrest are also studied numerically [Bj]. Active and living 



materials, such as self-propelled colloids and grains 
cells and bacteria 0, [l(| form a second class of nonequi- 
librium materials where glassy effects are reported and a 
dynamic arrest is observed at large density. Numerical 
studies of simple models of active matter in the dense 
regime have recently appeared (lll - H^ . similarities with 
the equilibrium glass problem were noted Q . More gen- 
erally, the issue of mapping the far-from-equilibrium dy- 
namics of active and self-propelled particles to an 'effec- 
tive' equilibrium problem is a recurring theoretical ques- 
tion [IM3, 

for which we are able to provide an ex- 
plicit solution in the particular context of the dynamics 
of driven glassy materials. 

Previous investigations either assume from the start, or 
demonstrate through detailed observations, that the ob- 
served dramatic change in dynamic properties is 'similar' 
to the glass transition of simple fluids at thermal equilib- 
rium [18J, in the sense that typical signatures of glassy 
dynamics (caging effect, two-step decay of time corre- 
lation functions, dynamic heterogeneity, etc.) are found. 
Theoretical studies exist for particular systems [5L [TlL fig| . 
and this broad set of experimental observations raises, we 



believe, a more general question that was not adressed in 
previous work. Can a glass transition occur far from 
thermal equilibrium, and how similar is it to the corre- 
sponding equilibrium phenomenon? 

Here we address this question from a theoretical view- 
point. We study the effect of non-thermal sources of en- 
ergy injection and dissipation on the slow dynamics of 
glassy materials. To obtain quantitative results, we fo- 
cus on a specific, yet representative class of glass models 
whose equilibrium behaviour can be studied analytically 
and is well understood [20j . following exactly the gen- 
eral scenario obtained in the context of random first or- 
der transition theory 2lj. Our main result is that glass 
transitions might exist even when driving forces are not 
thermal. Our results also indicate that the location of 
the transition continuously depends on the microsopic 
details of the dynamics, but that slow relaxation near 
dynamic arrest is insensitive to these details, thus sug- 
gesting that an 'effective' equilibrium glassy dynamics 
generically emerges at long times in nonequilibrium ma- 
terials close to dynamic arrest. 

To model the experimental situations described above 
we consider the dynamics of N degrees of freedom, x = 
{xi,i = 1---N}, representing for instance the position 
of grains or cells, interacting through the Hamiltonian 
%[x], which is supposed to display a glass transition at 
thermal equilibrium. The driven and active materials we 
wish to study share two important characteristics. First, 
they dissipate energy through internal degrees of freedom 
at a finite rate. Second, energy is continously supplied 
either by a global external forcing or by the particles 
themselves. To account for these effects, we study the 
following equation of motion: 



Xi{t)- 



ds 7d (t-s)x i (s) + 7 ^+r h (t) + f?(t) = 0, (1) 

UXi 



where we have included contributions from both a 
(white noise) equilibrium thermal bath satisfying the 
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fluctuation-dissipation theorem, (i]i(t) , r]j(s)) = 2TS(t — 
s)5ij and from nonequilibrium, colored driving and dissi- 
pative mechanisms represented by the active force /"(t) 
and the dissipation kernel jd(t), respectively. Equation 
(TTJ) is a standard theoretical model for the dynamics of 
active colloids and molecular motors far from equilib- 
rium [22i - [2^ | . It also represents a minimal model to an- 
alyze the physics studied in numerical treatments of ac- 
tive ll|, ll2J] and self-propelled particles 13[ , where parti- 
cles perform persistent random walks. It certainly misses 
some features of more complicated situations, such as 
complex alignement rules or particle anisotropy pjj . 

A more intriguing relation, valid when the friction (but 
not the noise) is Markovian, is with a Langevin pro- 
cess j2(| that mimics a quantum equilibrium bath, and 
reproduces some features of the true (operator-valued) 
quantum Langevin equation 27] ■ The parameters e a and 
T a used here are the counterparts of H and 

Equant 

We shall use this analogy below. 

As a first step we choose simple functional forms 
for the colored noise and dissipation terms. We use a 
Gaussian random forcing with mean zero and variance 
(/"(*)/"(«)> = 2F a (jt-8)6ij, where F a (t) = ^exp(-i), 
with r a the timescale of the slow forcing. For the dissi- 
pation we similarly choose 7d(f) = ^ exp(— which 
defines the timescale r^. With these definitions, ther- 
mal equilibrium is recovered either when = e Q = 0, or 
when the colored forces and friction satisfy the equilib- 
rium FDT [28], F a (t) = Tj d (t). We expect similar results 
for different functional forms for these correlators, as long 
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as they describe noises with finite correlation times. 

To make the problem tractable, we perform a mean- 
field approximation of the glass Hamiltonian. Our goal 
is to have a well-understood equilibrium starting point 
to isolate the influence of the non-thermal forces. We 
specialize our study to the spherical p-spin model, as a 
well-known representative microscopic model where the 
random first order transition scenario becomes exact [20j j . 
It is defined by 



u 



E 



(2) 



for continuous spins obeying the spherical constraint 
X], xj = N. In short, our strategy is to perform a mean- 
field approximation to the interactions in the equation of 
motion l|T]). while retaining realistic forms for the sources 
of injection and dissipation. As usual when dealing with 
glassy dynamics, our theoretical predictions strictly hold 
within the particular context of random first order tran- 
sition theory, but we expect them to have wider physical 
relevance, see Ref. [l8| for a broad theoretical overview. 
Another advantage of our approach is that it provides 
precise predictions which are then useful guides to com- 
puter simulations of more realistic models of active par- 
ticles. 

Because the Hamiltonian ([2]) is fully-connected, closed 
and exact equations of motion can be derived for the 
autocorrelation function C(t,s) — (xi(t)xi(s)), and for 
the autoresponse function R(t,s) = d{xi(t)) /drji(s): 



dC{t,s) 
ds 

dR(t, s) 
ds 

/*(*) 



H(t)C(t,s)+ dt'D(t,t')R(s,t') + dt'Z{t,t')C(t',s)+2TR(s,t), 



-fx(t)R(t,s)+ / dt'E{t,t')R(t',s) + 6(t- s), 

J —OO 

T+ I dt'[D(t,t')R(t,t') + Y,(t,t')C(t,t')], 

J —OO 

I 



(3) 



with the kernels D(t,s) = ^C p - 1 (t,s) + F a (t - s) and 
E(t,s) = E(p=llcP- 2 (t, s )R(t, s ) + 9 ^- s) . The last 
equation in Eq. enforces the spherical constraint. 

Technically, introducing colored friction and noise 
breaks detailed balance and introduces new physical 
timescales (t<2 and r a ) which compete and perturb the 
equilibrium dynamics of the system. This situation is 
reminiscent of the driven dynamics of the model stud- 
ied in previous work [29j j . where non- Hamiltonian driv- 
ing forces were introduced to model an applied shear flow. 
The crucial difference is the form of the driving terms, 
whose typical timescales in Ref. [2^] were that of the dy- 



namics itself, while here they relax with their own, fixed 
timescales Ta and r a (cf. Eq. ©). As a result, while the 
equilibrium glass transition was found to disappear in the 
presence of power dissipation of infinitesimally small am- 
plitude [29(, we show here that the glass transition may 
also survive the introduction of fluctuating forces, even of 
large amplitude, a situation which is closer to the effect 
of an oscillating field [3(| . 

The dynamic equations ([3]) are coupled integro- 
differential equations for which no general analytical so- 
lution can be obtained. Our strategy was to explore the 
region in parameter space where the dynamics becomes 
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FIG. 1: Slow dynamics with colored dissipation and noise, (a) Time correlators obtained from the numerical solution of Eqs. (|3]) 
for p = 3 at constant temperature: Increasing the amplitude of the dissipation ed with Td = 10 drives the system to the glass 
phase, (b) Decreasing temperatures for constant driving force with e a = 0.3, r a = 3. (c) Phase diagram of the system for active 
forces (T c < T^) and slow dissipation (T c > T^ q ), for various timescales Td,r a = 1,3, 10,30, 100 (from bottom to top). The 
vertical line indicates T,? q , dashed lines indicate the limits r<j,r a — ¥ 0, Eq. (|4j|. 



stationary, such that a driven steady state is reached and 
the system is in a 'fluid' phase with a finite relaxation 
time. We obtained numerical solutions of the dynamical 
equations, which motivated an analytic investigation of 
their asymptotic properties. 

We have adapted the numerical scheme of Ref. to 
obtain a solution of Eqs. ([3]) for a broad range of space 
parameters, using the following values for the time scales 
Td and r a : (1, 3, 10, 30, 100), in units of the bare friction 
coefficient which was set to unity in Eq. ([1]). We varied 
td and e a to fully explore the fluid part of the phase dia- 
gram. Numerically, we are limited by the stability of the 
numerical algorithm used to solve equations that do not 
have a formally causal form in the time differences [29| . 
We illustrate our results in Fig. [1] with the solution ob- 
tained in two limiting cases where either only dissipation 
or active driving forces are introduced. We do not find 
any added complexity when both terms are simultane- 
ously present with different amplitudes and timescales. 

A direct analytical solution only exists in the trivial 
limit where r^, r a — > 0, because Eq. ([T]) reduces to a stan- 
dard Langevin dynamics with white noise and memory- 
less friction. Equilibrium is then achieved at the rescaled 
temperature 



T 



T 



(4) 



Below, we argue that for finite correlation times of driv- 
ing and dissipation forces, a (more complicated) relation 
to the equilibrium situation only holds in the limit of 
very large time scales, a result which is far more sub- 
tle than the mapping in Eq. @. Equation (Q} implies 
in particular that response and correlation functions sat- 
isfy the fluctuation-dissipation theorem, TR(t) 
and that Eqs. © reduce to 



dCjt) 
dt 



dt w 2T J y ' dt' 



(5) 



which is mathematically equivalent [20( to the so-called 
F p _i schematic model derived in the context of the mode- 



coupling theory (MCT) of the glass transition [3l|. Its 



solution is known in great detail, and displays a dynamic 
singularity as the (rescaled) temperature is lowered to- 
wards the equilibrium value T° q . Near the dynamic tran- 
sition, the correlation function develops a two-step de- 
cay, with asymptotic time dependences that follow the 
behaviour described for discontinuous (or 'type B') tran- 
sitions within MCT [3l|. An expression for the critical 
temperature and scaling laws is obtained by performing 
a detailed mathematical analysis [31( of the equation de- 
rived by taking the long-time limit of Eq. ([5]) , 



1 



T 



fC s {t) + 7p / dt'D s {t-t') 



,.dC s (t') 1-q 



dt' 



T 



D s (t)=0, 
(6) 

where D s {t) = fCJT 1 ^)- 

Note that, in the general driven case with a nonequi- 
librium bath, one can only go from Eqs. ([3]) to Eq. §5§ 
by formally eliminating R(t — t') from the dynamic equa- 
tions. In this way, one only obtains a version of Eq. ([5]) 
with the memory kernel taking a form M [C] (t—t 1 ), which 
is a complicated nonlocal functional of C (t) (as opposed 
to the local functions like C p_1 (i)). 

Our central result is that the main features of this equi- 
librium glass transition robustly survive the introduction 
of a finite amount of non-thermal fluctuations driving 
the system far from thermal equilibrium, as illustrated 
in Fig. [TJ In both cases of colored noise or dissipation, 
we numerically find that time correlation functions dis- 
play a two-step decay reminiscent of the equilibrium be- 
haviour. This emerging glassy dynamics is characterized 
by a relaxation time that diverges upon approaching a 
dynamic transition, which we call a nonequilibrium glass 
transition. We find that the location of the transition 
is a continuous function of the driving mechanism. The 
data in Fig. [T] confirm the natural expectation that T c in- 
creases in the presence of the additional colored dissipa- 
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tion, T c (ed,Td) > T° q , while it decreases in the presence 
of a colored forcing, T c (e a , T a ) < T° q . The phase diagram 
determined numerically in Fig. [TJ; suggests a nearly lin- 
ear dependence of T c on either e a or eg, with a less pro- 
nounced dependence on the time scales (note that Td and 
Tc; vary over 2 orders of magnitude in Fig. [Tb). 

We can analytically rationalize the above findings, and 
obtain additional insight into the slow dynamics near 
nonequilibrium glass transitions, relying on the fact that 



the driving terms responsible for the explicit violation 
of detailed balance in the dynamical equations ([3]) have 
correlations that vanish at long times. Thus, a strong 
scale separation occurs when the structural relaxation is 
much larger than both Td and r , i.e. sufficiently close to 
the nonequilibrium transition. We seek an approximate 
equation of motion valid for stationary states in the limit 
of large times, r^, r Q <C t, corresponding to the approach 
and departure from the plateau: 



dCsjt) 
dt 

dR s {t) 



dt 



(/i - I s )C.(t) + I R D s {t) + I dt'V s (t - t')C s (t') + I dt' [D s (t + t')R s (t') + £.(t + t')C s (t')] , (7) 

(8) 
(9) 



-(/x-/ s ) J R s (f) + £ s (t)JH+ f dt'Z s (t-t')R s (t'), 

Jo 

POO 

T + il+ dt'[D s (t')R s (t') + E s (t')C s (t')), 
Jo 



r 



with £, 



?&^lCP- 2 R s , D s = fCT 1 , and the fol- 
lowing integrals Z s = / °° dfE/CO. I R = f™ dt'R f (t'), 
and Q = f™dt'[D f (t')R f (t') + Z f (t')C f (t% We de- 
fined the 'slow' functions C s (t) and R s (t) as the exact 
solutions of Eqs. (JTHSJ), while the 'fast' ones are defined 
by difference, e.g. Cf(t) = C(t) — C s (t), and decay over 
time scales that do not diverge at the transition. 

A crucial element of the dynamical equations ([7j|9|) gov- 
erning the long-time dynamics is that the terms respon- 
sible for the explicit breaking of detailed balance have 
disappeared. They appear very indirectly through time 
integrals over the short-time dynamical behaviour. As an 
immediate consequence, these equations can be consid- 
erably simplified because Eqs. ([7j [8]) reduce to the same 
equation if correlation and response satisfy 



Rs(t) 



Xs dCsjt) 

T dt ' 



(10) 



which defines the fluctuation-dissipation ratio, X s , or 
equivalently an effective temperature T e g = T/X s [32J. 
A similar ansatz holds in the long-time limit of the aging 
regime (33j, i.e. in the unperturbed glass phase, and in 
the equivalent limit of vanishing shear- like forces [29j . 
Combining Eqs. (ff] l9j fTCTjt . we obtain 



(t + |aw 



p 



f dt 'D s (t-t')^P--I R D s (t) 



0(11) 



where q represents the intermediate plateau height of 
C(t). This equation constitutes our main analytical 
achievement. It shows that, at sufficiently long times, 
the dynamical equation governing structural relaxation 



is equivalent to the one found for equilibrium relaxation, 
Eqs. ([5l [5]), showing that an equilibrium-like glassy dy- 
namics emerges out of nonthermal forces driving the dy- 
namics at short times. A similar conclusion was reached 
in Ref. @ for the specific case a driven granular fluid, but 
here we find that this form only holds once the fast mo- 
tion is averaged away. Even then, there remain, however, 
two important differences with the equilibrium case. 

First, the 'coupling' parameters determining the nu- 
merical value of the critical temperature are 'renormal- 
ized' by the microscopic details of the driving forces 
through time integrals over the short-time dynamics, as 
can be seen by directly comparing Eqs. ^ and (fTT|) . 
This explains our numerical finding that the location of 
the transition and the value of the plateau height con- 
tinuously depend on the details of the microscopic dy- 
namics, Fig. fX}:. Thus, an analytic determination of the 
locus of the nonequilibrium glass transition requires solv- 
ing not only Eq. (TTTT) , but also strongly nonuniversal fea- 
tures of the short-time dynamics. By contrast, because 
the long-time dynamics remains described by a discon- 
tinuous MCT transition, all universal features of the time 
correlation functions remain valid far from equilibrium, 
as found numerically in Fig. [TJ In particular, while the 
critical exponents describing the time dependence of cor- 
relation functions will depend on the driving forces, the 
well-known relations between them 31[ remain valid. 

Second, while Eq. (jlip simply involves the correlator 
C s (t) , as in equilibrium, the response function R s (t) does 
not obey the equilibrium fluctuation-dissipation relation, 
but only an 'effective' one, Eq. (fTO)) . Note that we pre- 
dict the existence of a nonequilibrium T e ff for slow de- 
grees of freedom even for the stationary fluid phase, not 
only deep into the glass as in Ref. [34j . This finding illus- 
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FIG. 2: Effective temperatures near nonequilibrium glass transitions, (a) Dissipative forces produce T c g < T (i.e. X a > 1). 
Parameters as in Fig. [Tk. (b) Active driving produces T e g > T (i.e. X s < 1). Parameters as in Fig. [lp. (c) Effective critical 
temperatures for parameters as in Fig. [Tb. 



trates that nonequilibrium glass transitions are conceptu- 
ally distinct from the equilibrium analog, but that a form 
of equilibrium-like glassy dynamics naturally emerges at 
long times. 

Our numerical analysis confirms the existence of effec- 
tive temperatures, see Fig. [5] As usual [lH, [33J, we rep- 
resent integrated response functions, x(i) = f*dt'R(t'), 
versus C(t) using t as a running parameter, which yields 
straight lines with slope —X s /T whenever Eq. (| 10[) holds. 
The straight segments corresponding to the first step 
of the relaxation have no reason, despite appearances, 
to be straight since the system is strongly driven. The 
data in Fig.[2]show that X s behaves differently if friction 
(X s (e d ,T d ) > 1), or forcing (X s (e a ,T a ) < 1) dominates 
the physics. The former represents an unusual situation 
where slow degrees of freedom appear to be colder than 
the bath [35| . Although Eq. ((U) cannot be used to predict 
the actual value of T e g in the general case, it correctly 
predicts its qualitative trends and thus provides a sim- 
ple physical argument for its variation with our control 
parameters. These two distinct cases are reminiscent of 
the distinction between adamant ('hot') and susceptible 
('cold') molecular motors of Ref. [19(. Note also the strik- 
ing similarity with the corresponding curves measured in 
active fluids [III EH ■ 

Finally, we summarize our numerical findings for X s 
measured at the dynamic transition into a 'phase dia- 
gram' (e,T* s = T c /X s ) in Fig. [2b. With friction alone, 
we confirm that T* s < T c (ed, Td) (i.e. X s > 1), but the re- 
sulting T* s is lower than the equilibrium value T^ q . The 
additional dissipation thus drives the transition temper- 
ature T c up, but it shifts the 'effective' critical temper- 
ature T* ff down. Active forces have the opposite effect, 
see Figs. EJxc. Somewhat disappointingly, the physical 
intuition [l9( that nonequilibrium transitions occur at 
shifted T c but constant T* s w T° q is not valid for our 
simple model: the effective temperature cannot be used 
to infer the nonequilibrium phase diagram in our model. 
This can also be seen analytically, as the mapping be- 
tween equilibrium and off-equilibrium dynamics is more 



complicated than simply replacing T by T c g, see Eqs. ([6l 

EETJ- 

Let us discuss the general picture of the (T, e) phase 
diagram, Fig. lc, in particular for the active forces. The 
red lines mark a limit between a liquid phase (higher T 
and e) and a glass phase where the dynamics performs 
aging. Our numerical method does not allow us to deter- 
mine whether this transition line continues up to T = 0, 
as the numerical solution of the equations becomes un- 
stable. One possibility is that the red lines continue and 
incide at the T = axis without qualitative changes, thus 
enclosing the glassy phase. An alternative is suggested 
by the close analogy with the quantum Langevin equa- 
tion described above. For the case of the (T, h) plane, it 
is well established (3(| that a continuous transition line 
(very much like our red lines) exists up to a point where 
the transition becomes thermodynamically first order all 
the way to T = (see Fig. 2 of Ref. j36[). For our model, 
it would mean that the transition crossed by decreasing 
e at constant (very) low temperature is first order: the 
system would vitrify suddenly from a liquid phase with- 
out any precursor increase in viscosity. Such an intrigu- 
ing phenomenon, seemingly consistent with the data in 
Ref. clearly deserves further study. 

Our study provides firm theoretical grounds to the 
emerging view that dynamic arrest in driven and ac- 
tive materials shares important similarities with the equi- 
librium glass transition, in particular regarding the be- 
haviour of time correlation functions and the emergence 
of an 'effective' thermal behaviour. Since our study was 
performed within a particular theoretical framework, it 
would be interesting to check the robustness of our find- 
ings beyond the realm of schematic glassy models, for 
instance using nonequilibrium mode-coupling types of 
approximations (37l . |38| . Further away from mean-field 
approximations, one could for instance also test in ac- 
tual measurements whether the role of 'activated pro- 
cesses' 18] remains the same in and out of equilibrium, 



and if the mode-coupling crossover is as relevant in driven 
systems as it is in equilibrium ones. 
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Regarding effective temperatures, a numerical study 
of a kicked granular fluid [5| reports a shift of the glass 
transition towards larger density with increasing activity, 
suggesting that drive dominates over friction in that case. 
Thus, we predict the existence of an effective tempera- 
ture T e s > T for this system. Such a result i s rep orted in 
numerical studies of active particle systems H^|, but 
these measurements were not performed close to dynamic 
arrest. Closer to our ideas are the experimental mea- 
surements of an effective temperature in an air-fluidizcd 
granular bed [39| . then used to infer a nonequilibrium 
equation of state. The reported deviations from equilib- 
rium behaviour (39| are highly reminiscent of our find- 
ings in Fig. [2j We believe that future studies of active 
particles at large density [l3j], or biologically driven sys- 
tems [23 | will result in more examples of nonequilibrium 
glass transitions as described in this work, where tests of 
our predictions about the universality of the slow dynam- 
ics and the emergence of effective temperatures could be 
performed. 
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